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Abstract 

Recent experiments on model catalysts have shown that Ag-Cu alloys have improved selectivity 
with respect to pure silver for ethylene epoxidation. In this paper we review our first-principles 
investigations on the (111) surface of this alloy and present new findings on other low index surfaces. 
We find that, for every surface orientation, the presence of oxygen leads to copper segregation to the 
surface. Considering the alloy to be in equilibrium with an oxygen atmosphere and accounting for 
the effect of temperature and pressure, we compute the surface free energy and study the stability 
of several surface structures. Investigating the dependence of the surface free energy on the surface 
composition, we construct the phase diagram of the alloy for every surface orientation. Around the 
temperature, pressure and composition of interest for practical applications, we find that a limited 
number of structures can be present, including a thin layer of copper oxide on top of the silver 
surface and copper- free structures. Different surface orientations show a very similar behavior and 
in particular a single layer with CuO stoichiometry, significantly distorted when compared to a 
layer of bulk CuO, has a wide range of stability for all orientations. Our results are consistent 
with, and help explain, recent experimental measurements. 

PACS numbers: 68.43.Bc, 81.65.Mq, 61.82.Bg 
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I. INTRODUCTION 



Processes such as heterogeneous catalysis and corrosion are dominated by chemical reac- 
tions that happen on the surfaces of materials These processes involve length and time 
scales that span several orders of magnitude: on one hand, molecular reactions occur on a 
sub-nanometer length and picosecond time scale, while, on the other hand, for surface phase 
transitions, catalysis, and corrosion the relevant time scale is of the order of microseconds 
or even seconds. Clearly this calls for a multiscale approach if one aims at describing, e.g., 
the full catalytic cycle or the full corrosion problem. 2j 

Our present study aims at obtaining a detailed knowledge of the atomic and molecular 
processes that underlie such phenomena. We are therefore interested in providing an atom- 
istic description of the transformations that take place at surfaces, in particular in the case 
of heterogeneous catalysis. The study of the formation or the breaking of a chemical bond, 
in any reaction, requires a quantum-mechanical description of the electrons. This calls for a 
first-principles description that does not rely on empirical parameters but rather just on the 
knowledge of the types of atoms present in the materials or molecules involved in the pro- 
cesses. It also requires the consideration of the quantum-mechanical kinetic energy operator 
in order to account for the shell structure of atoms, molecular orbitals, and chemical bonds. 
Among the electronic structure theories that have been developed to tackle such problems, 
density functional theory (DFT) j^, 4] has emerged as the most successful approach, given 
the very favorable ratio between accuracy and computational cost that this theory provides. 

As a first step toward modeling heterogeneous catalysis, an atomic description of the 
structure of the surfaces is required, For decades surface science techniques have been 
employed to study surface structures and transformations in ultra high vacuum (UHV) 
conditions and the knowledge thus obtained has been applied to "real life" catalytic processes 
that often take place at high temperature and high pressure regimes. While, in many cases, 
this approach has been extremely successful (see e.g. Ref. M), the pressure gap between 
surface science experiments and realistic applications inhibits, in several direct 
transfer: This is because the structure and composition of the surface can drastically change 
when going from UHV to realistic conditions Q, • A typical system that undergoes such 
a transformation is ruthenium: the catalytic activity toward CO oxidation under reactant 
pressures of the order of atmospheres has been in fact attributed to the formation of an 
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Ru02 oxide overlayer on Ru, rather than pure Ru (see Ref. [3] and references therein). 
Moreover, in this system the highest catalytic activity is reached in regions of the phase 
diagram corresponding to boundaries between different stable surface oxide structures, and 
the surface structure and stoichiometry at the steady state of catalysis is very different 
to any thermodynamically stable surface. Experiments on CO oxidation at Pt(llO) and 
Pd(lOO) conducted in high temperature and high pressure conditions have also shown 
how, under steady-state conditions, thin oxides form at the surface and they continuously 
evolve with time. These recent results have therefore highlighted both the importance of 
the environment and the non-equilibrium nature of catalysis: The catalyst surface cannot 
be viewed as a static object but rather as a system that, at the atomistic level, continuously 
evolves and fluctuates. 

It is therefore clear that the effects of temperature and pressure must be properly ac- 
counted for in order to obtain a reliable description of the surface structure and the processes 
happening on it. To this end, first-principles calculations, which combine accurate electronic 
structure methods with equilibrium thermodynamics, have been successfully employed to 
predict the stability of surface structures of several systems under high temperature and 
high pressure conditions, Moreover, such methods have been included in a multiscale ap- 
proach that employs kinetic Monte Carlo to simulate the dynamics of the full steady-state 
catalysis. In some cases the agreement with experiments can be quantitative, allowing an 
unprecedented level of insight into the concerted actions of a dynamical process such as 
heterogeneous catalysis. [10|] Hence, combining accurate electronic structure methods and 
non-equilibrium statistical mechanics techniques is now emerging as a possible way to tackle 
the multiscale nature of such processes. 

Furthermore, research carried out in the last few years on oxidation reactions on transi- 
tion metal catalysts have shown how the catalyst surface can present a number of possible 
structures, rather than a single low energy one. This suggests a scenario in which the cat- 
alyst dynamically evolves in time, fluctuating between structures with similar energetics. 
This picture has emerged, for example, in recent studies of CO oxidation at RuO2(110) [10,] 
and Pd(lOO) ll|: In both cases the highest catalytic activity is reached in regions of the 
phase diagram corresponding to boundaries between different stable surface oxide struc- 
tures. Experimentally, it has been found through scanning tunneling microscopy (STM) 
measurements conducted under realistic catalytic conditions for CO oxidation at Ft (110) 
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and Pd(lOO) [8| that, under steady-state catalysis, the thin oxides present at the surface 
continuously evolve with time. These recent results have therefore stressed how the cata- 
lyst surface cannot be viewed as a static object but rather as a "living" system that, at the 
atomistic level, continuously evolves due to various processes such as adsorption, desorption, 
association, dissociation, diffusion and to the competition between structures with similar 
energetics. 

The use of bimetallic catalysts has been the focus of much work in the field of heteroge- 
neous catalysis, |l2| since the catalytic activity and selectivity of a metal can be modified 
substantially by alloying with another metal. For example, geometrical and electronic ef- 
fects obtained by varying the alloy composition may play a crucial role in determining the 



properties of the catalyst. 



13| One appealing aspect of bimetallic alloys is the possibility of 



rationally designing the catalytic properties of the material by changing the alloy elements 
and composition. To this end, first-principles calculations have been shown to be a poten- 
tially useful tool for screening among pools of possible alloys and for extracting trends and 
therefore gaining insights into the functioning of the alloy. isj. In view of the above 
discussion, though, the effects of high pressure and temperature and the role of the dynam- 
ical evolution of the catalyst must be carefully accounted for, in order for the theoretical 
simulations to have predictive power. 

A catalytic process in which the use of bimetallic catalysts could provide a significant 
benefit is ethylene epoxidation, which is one of the most important processes in the chemical 
industry. Here, an ethylene molecule (C2H4) is oxidized to ethylene oxide (C2H4O), which is 
used in the synthesis of ethylene glycol, polyester fabrics, surfactants, and detergents. The 
catalyst for this reaction must promote the formation of ethylene oxide and at the same 
time inhibit the competing reaction, total oxidation to CO2, which is thermo dynamically 
more favourable. The selectivity toward the formation of ethylene oxide is therefore a key 
factor. In industrial applications, the catalyst used for this important chemical reaction is 
silver, and the reaction is carried out at atmospheric pressure and temperatures T=500-600 
K. [16| Both the selective (formation of ethylene oxide) and unselective (total oxidation to 
CO2, which occurs via the formation of acetaldehyde) paths have been found to originate 
from a common intermediate, an oxametallacycle, in which the molecule is bound to both a 
chemisorbed oxygen can react to form 

either ethylene oxide or acetaldehyde with similar activation barriers. On the basis of both 
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first-principles calculations [l9j and experiments [20] it has been suggested that if an Ag-Cu 
alloy, rather than pure Ag, is used as a catalyst, the selectivity toward ethylene oxide will 
be improved. 

In several recent works the mechanism behind the enhanced selectivity due to alloying 
with copper has been addressed. [21I, 22] The composition and structure of the catalyst 
surface, though, has not been thoroughly investigated. Given the discussion presented above, 
it is clear that the mechanism of this catalytic reaction can be strongly influenced by the 
atomistic details of the surface. In particular, given the conditions of operation of the Ag-Cu 
catalyst in realistic applications, the formation of surface oxides cannot be ruled out. It has 
been established that silver can form several surface oxide structures 23|, |2J, |25| and from 
the phase diagram of the Cu-0 system, considering the temperature and oxygen pressure 
used in practical applications, but neglecting the effect of ethylene, the formation of CuO is 
expected, therefore suggesting that copper might also oxidate. 

In this paper we review our work on the first-principle investigations of the structure of 
the Ag-Cu(lll) surface under varying oxygen pressure and temperature conditions, and we 
present new results on other low-index surfaces. In our calculations we consider the alloy 
surface to be in thermodynamic equilibrium with an atmosphere of pure oxygen. The effect 
of the presence of other reactants such as ethylene is not investigated in this work. Under 
the conditions of steady-state catalysis it is therefore likely that the surface of the catalyst 
might be modified with respect to the stable structures found in a pure oxygen atmosphere. 
Hence our work can be thought of as a first step in gaining insight into the possible structures 
of the alloy, an essential prerequisite for modeling the full catalytic process. The paper is 
organized as follows: In Sec. [TTlwe briefly review the theoretical background relevant for this 
work, present the definitions of the quantities that will be used throughout the paper and 
discuss the approximations and assumptions that underpin our methodology. In Sec. [1111 we 
report our results and discuss their relevance. Section HVl summarizes the main findings. 

II. CALCULATION METHOD 

The density functional theory (DFT) calculations presented in this work are performed 
using the generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) {2^ 
for the exchange and correlation functional. The electron-ion interactions are computed us- 
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ing ultrasoft pseudopotentials, [27|, l28| including scalar relativistic effects, and the energy 



cutoff for the plane wave expansion is 27 Ry (200 Ry for the charge density cutoff). The 
Brillouin zone is sampled using a Monkhorst-Pack mesh, 2^ broadening the Fermi surface 



according to the Marzari-Vanderbilt cold-smearing technique, [30| using a smearing param- 
eter of 0.03 Ry. A 12x12x1 k-point mesh is used for the (1 x 1) surface unit cell of the 
(111) and (100) surfaces, while a 12x9x1 mesh is used for the (1 x 1) surface unit cell of the 
(110) surface. For larger unit cells, the k-point mesh is scaled accordingly to give us iden- 
tical k-point sampling. All the calculations are performed using the PWscf code contained 
in the Quantum-ESPRESSO package, jsi] We use a slab geometry to model the surfaces, 
using 4 Ag atom layers for the (111) surface, 6 for the (100) and 6 for the (110), on top of 
which the adsorption structures are created, on one side on the slab only. For large unit 
cells, in particular for the (4 x 4) structures on the (111) surface, we reduced the number 
of Ag layers to three. The bottom two layers are fixed to the bulk positions, while all other 
atoms are allowed to relax until the forces are less than 0.001 Ry/Bohr (0.025 eV/A). A 12 
A vacuum layer is used, which is found to be sufficient to ensure negligible coupling between 
periodic replicas of the slab. Further details of the calculations presented in this work can 



be found in one of our previous publications 



32|. 



To help the discussion presented in the following section, we now define several quantities. 
Considering the adsorption of atomic oxygen on a clean Ag surface, we define the average 
binding energy per oxygen atom as 

No 

where A^o is the number of oxygen atoms in the unit cell, E'^^^^, ^siab j^o ^^^q 
total energies of the adsorption system, the Ag slab and half the oxygen molecule, i.e. 
E'^ = 1/2Eq^^^. With this definition, a positive binding energy indicates that the adsorption 
of an oxygen atom is exothermic with respect to oxygen in molecular gas phase form. Here 
we point out that the significant error (~ 0.54 eV) in the oxygen molecule binding energy 
introduced by the GGA-PBE approximation (l/2£'^^ = 3.10 eV, while the experimental 
value is 2.56 eV ^J) introduces a large error bar in the determination of the oxygen chemical 
potential. This error, however, is partially be compensated by the analogous GGA-PBE error 
in the description of oxygen chemisorbed on the metal surfaces. 

To account for the effects of temperature (T) and pressure (p) we employ "a6 initio 
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fly, 



atomistic thermodynamics" , |5|,l3j, l35l . l36l . l37l | which allows the determmation of the lowest 



energy structures as a function of T and p. The surface is considered to be in contact with 
an oxygen atmosphere that acts as a reservoir, therefore exchanging oxygen atoms with the 
surface without changing its temperature and pressure (i.e. its chemical potential). 
The change in Gibbs free energy is calculated as 

AG(/io) = ^(G^/^s _ ^siab _ A^r^^^^^ _ n^^^^ ^ (2) 

where G^/^^ jg \^]^q f^-gg energy of the adsorbate/substrate structure, G^^^^ is the free energy of 
the clean Ag slab, AA^Ag is the difference in the number of Ag atoms between the adsorption 
system and the clean Ag slab, /lAg and /io are the atom chemical potentials of Ag and O. 
The change in Gibbs free energy is normalized by the surface area A to allow comparisons 
between structures with different unit cells; we will refer to this quantity as the "change 
in Gibbs surface free energy of adsorption" , or simply "surface free energy" . The chemical 
potential of Ag is taken to be that of an Ag atom in bulk Ag, therefore assuming that the 
slab is in equilibrium with the bulk, that acts as the silver reservoir. The oxygen chemical 



potential depends on temperature and pressure according to _38| 

/^o(T,p) = i[Eg'fi(T,p°)+/io.(T,/) + A;BT/n(^)] . (3) 

Here p° is the standard pressure and jlo^{T,p^) is the chemical potential at the standard 

391 (the choice made in 



pressure, which can obtained either from thermochemical tables 
this study) or directly computed. 

While the free energy should include terms that depend on the vibrational modes of the 
system, as well as configurational disorder, in our treatment we neglect such contributions. 
Especially when comparing systems with different stoichiometries, such terms can be non- 
negligible. For our system, however, through a simple estimate based on the Einstein model 
and approximating the phonon density of states by just one characteristic frequency for 
oxygen, we have shown [32] that neglecting the vibrational contributions to the free energy 
does not alter the main findings of our work. Also the configurational entropy has been 



found, through a simple estimate, not to contribute significantly to the free energy. [32 1 
Hence, within this model, the only term that depends on T and p is the oxygen chemical 
potential. As a consequence, the free energies G^^^^ and G^^'^^ are identified as the total 
energies E^/^s and E'^'^'^. 



If we now consider the presence of copper impurities in the silver surface, we need to 
modify the definition of surface free energy. As we will show later in the paper, we deal 
with structures in which an overlayer containing O, Ag and Cu is adsorbed on a clean Ag 
surface. If we consider Cu to be in equilibrium with a bulk Cu reservoir, the definition of 
surface free energy becomes: 

AG(/io) = - G^'^^ - AA^Ag/iAg - Ncuf^cn - iVo/^o) , (4) 

where G"-'/^"/'^^ is the free energy of the adsorption system, A^^cu is the number of Cu atoms 
and ficu the copper chemical potential. Since in this work, as we will show in Sec. IIIIl we 
are interested in identifying the structures belonging to the convex hull of the free energy 
vs. copper content curve, which depends on the curvature of such curve, rather than the 
absolute surface free energy, the choice of the Cu chemical potential is arbitrary, since the 
surface free energy depends linearly on it (see Eq. H]). We now define the average oxygen 
binding energy in the whole structure as 

^O/Cu/Ag ^ _ 1 [^O/Cu/Ag _ (^slab ^ AA^Ag/iAg + A^c./iCu + NoE'')] . (5) 
iVo 

Using the above-mentioned approximations, we can express the change in Gibbs surface free 
energy in terms of the average oxygen binding energy: 

AG(A/xo)= -^(AToE^/^^/^^ + AToA/xo) , (6) 

where the oxygen chemical potential is now referenced with respect to half an isolated O2 
molecule: A/io = /^o - \Eq^^^. 

In writing the change in free energy as in Eq. ([6]) we are neglegting the contribution to the 



free energy coming from the configurational entropy. It has been shown 32] that indeed this 
approximation introduces negligible effects for the range of temperatures we are interested 
in. 

III. RESULTS 

A. Oxygen-induced segregation of copper to the surface 



Due to their significant (13 %) lattice mismatch, copper and silver are known to be almost 
completely immiscible in the bulk at temperatures of interest for industrial applications. In 
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several cases, however, elements which are immiscible in the bulk have been found to form 
stable two-dimensional alloys at the surface 



401 . These include Na and K deposited on 

' ' L n n 

Al(lll) and (001), [41] Au on Ni(llO), |42|] Ag on Pt(lll) [43^ and Sb on Ag(lll). j44| 
Indeed, the presence of such a two-dimensional alloy has been invoked also for the Ag-Cu 
system in the work by Linic et al. [19|] 

Considering the (111) surfaces of Ag and Cu, due to the fact the surface energy of silver is 
considerably lower than the one of copper (experimental values: 7Ag(i") = 0.078 eV/A^ |45| 
and 7Cu(iii) = 0.114 eV/A^ |46|), it is unlikely for copper impurities in silver to be exposed 
on the surface. We verified this by computing the total energy of a (111) silver slab in 
which every fourth Ag surface atom is substituted with a Cu atom. By changing the layer of 
substitution, we found that the most favorable position for a Cu impurity is one layer below 
the surface. We also found that an increase in the copper content on the surface leads to 
an increase in surface energy. This is true for all the three surface orientations considered 
here, i.e. (Ill), (100) and (110). The results are shown in Fig. [T]as the curves with black 
dots. The positive slope of the curves indicate that it is unfavorable for Cu to migrate to 
the surface. The black curves in Fig. [1] have an almost negligible convexity, indicating a very 
small tendency to form a surface alloy. The configurational entropy contribution to the free 
energy, neglected in the present approach, would further enhance the stability of an alloy 
at finite temperature. At T=600 K we estimate this term to contribute 5 meV/A^ to the 
surface free energy. 32| We therefore find that for this system, when Cu is confined to the 
first layer and in the absence of oxygen, there is a weak tendency to form a two-dimensional 
alloy. Similar results have been found for an Ag-Cu alloy on the Cu(lOO) surface. 

Cu thin-film growth on Ag(lll) in ultra high vacuum (UHV) has been recently studied 



experimentally. 
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491] It was found that at room temperature, upon deposition of Cu on 



Ag(lll), Cu forms islands that are encapsulated by one monolayer of Ag. These findings 
agree with our theoretical results which show that the most favourable configuration for Cu 
is to sit one layer below the Ag surface layer. 

When oxygen is introduced into the system, the picture changes completely: We find that 
the presence of oxygen chemisorbed on the alloy surface has the remarkable effect of reversing 
the slope of the (black dots) curves for the surface energy versus Cu surface composition, 
shown in Fig. [TJ For these calculations we consider an Ag-Cu alloy in the first layer of the 
Ag slab, with oxygen adsorbed in the most favorable adsorption sites. By fully relaxing 
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FIG. 1: (Color online) Change in Gibbs surface free energy at T = K (cf. Eq. [6|) as a function of 
Cu concentration in the first layer of the surface. The oxygen coverages considered are indicated 
in the legend in monolayers (ML). The energy zero corresponds the surface energy of the clean 
Ag(lll) surface. The top left panel refers to the (111) surface, the top right to the (100) surface 
and the bottom one to the (110) surface. 

the structure, we calculate the change in surface free energy at T = K as a function of 
the Cu content in the first layer and the oxygen coverage. The results are shown in Fig. [1] 
for the three surfaces considered in this work. We can see that for the (111) surface the 
presence of at least a quarter of a monolayer (ML) of oxygen induces Cu to segregate to 
the surface, i.e. it changes the slope of the curve from positive (for an oxygen coverage 
relative to the underlying Ag(lll) lattice 9o < 0.25 ML) to negative (for Oq ^ 0.25 ML). 
The driving mechanism here is the strong affinity between oxygen and copper, which more 
than compensates the unfavorable surface energy of Cu with respect to Ag: The adsorption 
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energy (cf. Eq.([T])) of oxygen at 0.25 ML coverage on Ag(lll) and Cu(lll) is 0.38 eV/atom 
and 1.57 eV/atom, respectively. We also see that in the absence of copper (i.e. the left end 
of the graphs in Fig. [T]), the most stable structure is the one with 6*0 = 0.25 ML, i.e. the 
lowest O coverage. At higher oxygen coverages the oxygen-oxygen repulsion overcomes the 
formation of 0-Ag bonds, leading to an increase in surface energy with the oxygen coverage. 
At the opposite end of the graph, where a full monolayer of copper is present in the first 
layer of the slab, the surface energy decreases with the oxygen coverage: This is due to the 
contribution to the surface energy of the formation of strong 0-Cu bonds. The (100) and 
(110) surfaces display a behavior qualitatively similar to the (111) surface, but the transition 
to a negative slope takes place at a higher oxygen coverage, 0.50 ML and 0.75 ML for the 
(100) and (110) surface respectively. 

A study similar to the one shown here has been recently presented for the AgsPd alloy. [5C 
Also in that system the presence of oxygen has the effect of inverting the segregation profile: 
While oxygen poor conditions favour Ag-terminated surfaces, in oxygen rich conditions Pd is 
driven to the surface, due to the strength of the Pd-0 bond compared to the Ag-0 one. The 
Pt-Ru alloy displays analogous properties as well, 5l| where the more noble Pt is exposed 
in absence of oxygen, while the strong Ru-0 bond leads to Ru segregation at the surface 
in oxygen rich conditions. All these systems therefore display the same trend: In absence 
of oxygen the element with the lowest surface energy is exposed, while in the presence of 
oxygen the element with the higher affinity with oxygen segregates at the surface. 



B. Stable structures on the Ag-Cu alloy surfaces 

Having established that in the presence of oxygen, for all the three surface orientations 
considered, there is a clear tendency for Cu to segregate to the surface, in the following 
we assume that Cu is confined to the surface and we investigate what are the most stable 
surface structures as a function of temperature and pressure. Given the fact that copper 
might oxidize under the typical conditions of oxygen pressure and temperature in which the 
catalyst is supposed to operate, we also consider surface structures that are derived from 
the corresponding bulk copper oxide geometries. Here we point out that the real catalyst 
operates in the presence not only of oxygen but also of ethylene, which acts as the reducing 
agent. Since in our calculations we do not consider the presence of ethylene, the stability of 
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FIG. 2: (Color online) Top and side view of some of the low energy structures considered in this 
work: (a) (2x1) added row structure on the (110) surface, (b) "p2" on the (111) surface, (c) 
Ag/0(2x2) on the (100) surface, (d) 1 layer of CuO on the (111) surface, (e) 1 layer of CuO on 
the (100) surface, (f) 1 layer of CuO on the (110) surface. The red (black) atoms represent oxygen, 
the blue (dark grey) ones represent copper and the grey (light grey) represent silver. 

the structures we study here could be altered by accounting for the effects all the reagents. 

We study, for each of the three surface orientations, three types of structures: (i) 
chemisorbed oxygen on the Ag-Cu alloy in the first layer of the Ag slab (as described in the 
previous section), (m) structures derived from copper (I) oxide CU2O, whose structure can be 
visualized as trilayers of 0-Cu-O piled up on top of each other and {iii) structures derived 
from copper(II) oxide CuO. 

For the first set of structures we employ (2x2) unit cells and label the structures using 
the notation Oa;ML/Cu?/ML, where x and y are the content of O and Cu, expressed in 
monolayers with respect to the underlying Ag surface. In the case of the (100) and (110) 
orientations, we also study the missing row and the added row reconstructions, respectively, 
allowing for Cu alloying in the first layer (see Fig.[2](a)). For the (100) surface, in the absence 



12 



of Cu, also the c 
reported in Ref. 



'4x6) reconstruction has been considered, using the data for the energetics 
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For the second set of structures, the geometries are derived from bulk CU2O. We use 
various film thicknesses, from 1 up to 5 layers. For each surface orientation, forcing the 
CU2O film to be commensurate to the underlying Ag surface , results in a 7% compression 
of the CU2O film area compared to the bulk value. In the case of the (111) orientation, we 
use the label "p2" or "p4" depending on whether the periodicity of the structure is (2x2) 
or (4x4) with respect to the clean Ag surface. Of particular interest are the "p2" and the 
"p4-OCu3" structures. The atomic geometry of the "p2" structure is shown in Fig. [2](b). 
Here each Cu is linearly bonded with two O atoms and each O atom is bonded to three 
Cu atoms, in a ring-like structure similar to the ones proposed for thin oxide-like layers on 
Ag(lll)[23] and Cu(lll) 53|]. The "p4-OCu3" is a p4 structure in which an OCus unit has 
been removed. For the (100) orientation we use a {y/2x ^/2) unit cell and for the (110) 
surface we use a (2x2) cell. None of these structures on the (100) and (110) surfaces, as we 
will show later, are found to be present in the phase diagram. 

For the third set of structures, we consider thin layers of CuO-like configurations, with 
a thickness up to 5 layers. In this case we must bear in mind that bulk CuO is poorly 
described with DFT-PBE: CuO is a strongly correlated antiferromagnetic semiconductor, 
with a monoclinic structure (a=4.65 A, 6=3.41 A, c=5.11 A, (3=99.5°), where Cu is linearly 
bonded to 2 O atoms and O is bound to 2 Cu atoms. 5J] DFT-PBE, on the other hand, 
predicts CuO to be a metal with an almost orthorhombic structure (a=4.34 A, 6=4.01 A, 
c=5.22 A, /3=92.2°), in which each O atom is tetrahedrally coordinated to 4 Cu atoms and 
each Cu atom is bonded in a square planar geometry to 4 O atoms. The predictions for 
thick films of CuO must therefore be regarded as qualitative. The computed formation 
energy per oxygen atom in bulk CuO is 1.23 eV (experimental value 1.63 eV [55|) and the 
Cu-0 bond length is about 1.97 A (experimental value 1.67 A). Going beyond the DFT-PBE 

culations for CuO and CU2O using GGA+U plus GqWo, [56| 



57! ] are in progress. Given the fact that CuO does not have 



description is a major task. Ca^ 
as well as other methodology, 
a cubic structure, forcing the CuO thin layers to be commensurate with the underlying Ag 
substrate leads to distortions that depend on the orientation of the surface. In our study 
we consider (2x2) cells for the (111) and (100) surfaces and (2x1) cells for the (110). This 
leads to a compression of the Cu-0 bonds that, for every orientation, is of the order of 3%, 
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and the relaxed geometries of the structures strongly depend on the surface orientation (see 
Fig. [2](d-f)). We find that one layer of CuO is particularly stable for every orientation. For 
this structure we calculate the formation energy per oxygen atom to be 1.16 eV for the (111) 
surface and 1.26 eV for both the (100) and (110) surfaces. 

In addition to the structures considered in this work, it is likely that other structures 



involving thin oxide- 
works on 0/Ag(lll) 



ike layers with similar surface energies exist for this system. Recent 



23 



m 



25l | and 0/Pd(lll) [58] have shown that for these systems a 



multitude of structures with similar energetics can be found. 



C. Thermodynamic diagrams for the Cu-O/Ag system 

In Fig. [3] we report the change in Gibbs surface free energy of the structures considered as 
a function of the change in oxygen chemical potential. The three graphs correspond to the 
three surface orientations considered in this work. We also show at the bottom of the plot, 
for three values of temperature (300, 600 and 900 K), the value of pressure corresponding 
to the chemical potential shown in the abscissa of the plot. 

As we can see from Eq. (|6]), the slope of the lines in Fig. [3] is proportional to the oxygen 
coverage, i.e. the higher the oxygen content the steeper the line. In all the three graphs, 
the vertical dashed line at A/xq = —1-26 eV represents the change in chemical potential 
above which Cu oxidizes to CU2O, which, as one can show, corresponds to the computed 
heat of formation of CU2O. The vertical dashed line at A/iQ = —1.23 eV, on the other hand, 
represents the change in chemical potential above which Cu oxidizes to CuO. 

In the case of the (111) surface, we can see that in the range of values of the oxygen 
chemical potential considered (—1.6 < A/xq < —0.4 eV) the thermodynamically most stable 
structures are: Pure Ag for A/xq < -1-43 eV, the surface oxides "j92" and "p4-OCu3" 
(almost degenerate) between —1.43 < A/io < —1.26 eV, bulk copper(I) oxide CU2O for 
-1.26 < A/io < -1-23 eV and bulk copper(II) oxide CuO for A/io > -1-23 eV. 

For the (100) surface we find that no surface oxides are stable. For A/iq < —1.26 eV 
pure Ag is the stable structure, while above that value of oxygen chemical potential, as in 
the (111) surface, we find the two bulk copper oxides to be stable. 

For the (110) surface, on the other hand, we find that in the range —1.46 < A/io < -1-26 
eV a three layer thick CU2O surface oxide is the stable structure, while at lower oxygen 
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FIG. 3: Thermodynamic diagrams (i.e. change in surface free energy versus oxygen chemical 
potential) for the Ag/Cu/0 system. The top left panel correspond to the (111) surface, the top 
right one to the (100) surface and the bottom one to the (110) surface. The vertical lines separate 
the regions where different structures are stable. 



chemical potential the stable structure is pure Ag and at higher oxygen chemical potential 
the bulk copper oxides become stable. 

It is interesting to note that none of the alloyed chemisorbed structures (the structures 
labelled as OxML/Cu|/ML) are thermodynamically stable for any of the three surfaces, 
while there are regions in which two-dimensional surface oxides are stable, both on the 
(111) and the (110) surfaces. This situation is similar to what has been found for the 
0/Cu(lll)system, |53|] where it was argued that copper oxidation does not proceed via 
ordered chemisorbed structures, at variance with other transition metal surfaces. We also 
find that none of the 0/Ag structures, i.e. structures containing no Cu (not shown in Fig. [3]), 
are thermodynamically stable in the range considered here. This is not unexpected, given 
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the weaker strength of the Ag-0 bond with respect to the Cu-0 one. 



D. Surface phase diagrams 

Having identified, among those considered here, the thermodynamically most stable or- 
dered structures for the infinite (111), (100) and (110) surfaces, we now consider the situation 
in which we have a well defined Cu surface content in each of those surfaces. For a finite 
system (e.g. an Ag particle) of known dimensions and well defined Cu content, if we assume 
that all the Cu present in a spherical nanoparticle segregates and homogeneously distributes 
on the surfaces of the particle in a site of the Ag fee lattice, given the size of particle and the 
total Cu content, we can estimate the Cu surface content. For example in a 50 nm spherical 
particle with a 1 % Cu content (typical values of particle size and copper content of Ag-Cu 



alloys used in experiments in which the effect of Cu impurities was investigated j20|) the 
(average) surface Cu content is 0.35 ML, where 1 ML corresponds to the first layer of the 
particle being pure copper. Experiments carried out on Ag-Cu alloys with 0.1-1% Cu con- 
tent at ~ 500 K have shown that in an oxidizing environment the Cu surface content is in 



the range 0.1-0.75 ML. 



20| 



In the following, we will consider the three surfaces investigated in this work and use 
the results shown in Fig. [3] to predict which structures will be present on the surfaces as a 
function of the copper content and the ox yge n chemical potential. We also exploit the vast 



literature available for the 0/Ag system [23|, |2J] to include the most stable structures for 
the system in the absence of copper. 

As an example, we show explicitly the case for A/iq = —0.80 eV, corresponding to an 
oxygen pressure of 10~'^ atm. at the temperature of 600 K. Figure H] shows the change in 
surface free energy at A/io = —0.80 eV of all the structures we have considered on the (111) 
surface (a total of 55 structures, including the 0/Ag ones) as a function of the Cu surface 
content. By constructing the convex hull, i.e. the curve obtained by joining those structures 
that are most stable with respect to linear combinations of structures at other compositions 
that would yield the same total composition, we identify the structures (indicated by the 
red dots in Fig. Hj) that are stable against phase separation into any two other structures. 
Although in Fig. |3]we show only the portion of the diagram up to a Cu content of 1 ML, 
we have computed structures including up to five CuO and CU2O layers in order to be able 
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FIG. 4: (color online) Construction of the convex hull (indicated by the line) for the case of 
A^o = —0.80 eV. Only the portion up to 1 ML is shown in the figure. The energy zero corresponds 
to the clean Ag(lll) surface. 

to extrapolate the behavior of bulk CuO. 

For a specific Cu content, e.g. 0.25 ML, the convex hull plot helps to predict that a 
mixture of pure Ag and "p4-Cu30" (rather than a single ordered structure with 0.25 ML 
Cu content) will be present on the Ag(lll) surface at A/iQ = —0.80 eV, in a ratio given by 
the lever rule, i.e. by the mass conservation law. By repeating this scheme for all the values 
of the oxygen chemical potential in the range considered in this work, we build a phase 
diagram as a function of the oxygen chemical potential and the copper surface content. This 
is shown in Fig. [5] for the three surface orientations considered in this work. We find that 
the convex hull, for high Cu content, always includes the structure with the largest number 
of CuO layers (or CU2O layers at lower oxygen chemical potential), i.e. bulk oxide formation 
is favored beyond a certain Cu surface content. We therefore cut the plot in Fig. [Slat a Cu 
composition of 1.50 ML, where it is understood that the rightmost structure is "bulk" oxide. 
However, care must be taken in the meaning given to such "bulk" structure in the context 
of a finite system as the one we are modelling here. We interpret the numerical evidence 
of the presence of the bulk structure in the convex hull as a tendency for the formation of 
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thick patches of either CuO or CU2O. 

In the phase diagram shown in Fig. [5] we have exphcitly written the two structures 
coexisting in a number of regions around the values of interest of oxygen chemical potential 
and copper content. Around A/iq = —0.60 eV, for all three surface orientations, we find that 
the stable structures include clean Ag (on the (111) surface) or Ag-0 structures (Ag/0(c4x6) 
on the (100) surface and Ag/O*^^ on the (110) surface) at low Cu content, CuO(lL) for 
intermediate levels and bulk CuO at high Cu content. The fact that these features are 
general and do not depend strongly on the surface orientation allows us to draw some general 
conclusions. First of all, a 2D-alloy on the catalyst surface is never stable, for any of the low 
index surfaces examined here. The strength of the Cu-0 bond is the factor determining the 
absence of any surface alloy and the presence, on the other hand, of thin or thick layers of 
copper oxide. Secondly, a thin layer of CuO is particularly favourable for all orientations, 
making this structure a possible candidate to explain the high selectivity of the Ag-Cu alloy 
for ethylene epoxidation. 

The picture emerging from our results is consistent with recent experiments performed on 
the Ag-Cu system under catalytic conditions. [59|] In these experiments Ag-Cu nanopowders 
(~ 100 nm in diameter and 2.5% Cu) were used as the catalyst for ethylene epoxidation at 
T = 520 K and p = 0.5 mbar, corresponding to A^o = —0.69 eV. Through a combination 
of in-situ XPS, and Near Edge X-ray Absorption Fine Structure (NEXAFS) measurements, 
thin layers of CuO and CU2O are found be to present on the surface, while no signs on a 
surface alloy are detected. Areas of clean Ag are also present on the surface, in agreement 
with our theoretical calculations. 

These results also show that the simple 2D-alloy structure adopted in Ref. 



19 



for the 



Ag-Cu surface to model theoretically the ethylene epoxidation reaction is not stable and is 
significantly different from what we expect to be relevant under high pressure conditions. 



The model used in Ref. 



19l assumes a (2x2) periodicity for the Ag(lll) surface in which 



one atom out of four is replaced with copper, and oxygen is chemisorbed on it. Our results 
suggest on the other hand that, for this surface, a model that includes, depending on the Cu 
surface content, clean Ag(lll), patches of oxide-like thin layers and thick layers of CuO, is 
a more appropriate model for the Ag-Cu catalyst (111) surface. Moreover, we have shown 
that this picture is valid also for the other low index surfaces. 
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FIG. 5: Surface phase diagrams showing the structures belonging to the convex hull as a function 
of the Cu surface content and thechange in oxygen chemical potential. The top left panel refers to 
the (111) surface, the top right to the (100) surface and the bottom one to the (110) surface. 

IV. CONCLUSIONS 



By means of ab-initio atomistic thermodynamics we have investigated the low index sur- 
faces of the Ag-Cu alloy in equilibrium with an oxygen atmosphere, accounting for the 
effects of temperature and pressure. For all the surface orientations considered in this work, 
the presence of oxygen on the surface leads to copper segregation, due to the strength of 
the 0-Cu bond relative to the 0-Ag one. We have studied, for every orientation, surface 
structures that include chemisorbed oxygen on a 2D Ag-Cu alloy as well as CU2O and CuO 
layers, both in the single layer limit and in the bulk limit. Through the construction of 
the oxygen chemical potential dependent convex hull, we have predicted, as a function of 
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the surface copper content, temperature and pressure, the combinations of structures that 
are most stable. The phase diagrams constructed this way show remarkable similarities for 
different surface orientations. For the (111) surface, around the temperature and pressure 
of interest for practical applications and depending on the copper surface content, our re- 
sults indicate that clean Ag(lll) and thin copper oxide-like structures (p4-Cu30, p2 and 
CuO(lL)) can coexist. For the (100) orientation the Ag/0(2x2) and the CuO(lL) struc- 
tures are predicted to be present on the surface, while on the (110) surface the Ag/O™"^ 
and the CuO(lL) structures are present in the phase diagram around the region of interest. 
Therefore a combination of Ag/ O structures and thin copper oxide layers are predicted for 
all surface orientations. A single layer of CuO, in particular, has a wide range of stability 
on all surfaces. The model for the surfaces of the Ag-Cu catalyst that emerges from our 
calculations differs substantially from the one proposed in earlier works aimed at investigat- 
ing the effect of copper impurities in silver on the mechanism of ethylene epoxidation. Our 
results suggest that to gain insight into the full catalytic cycle one should consider oxygen 
species such as those in the proposed structures, rather than just oxygen chemisorbed on 
the clean alloy surface. 

Here we have limited ourselves to the study of the thermodynamic equilibrium between 
the Ag-Cu alloy and an oxygen atmosphere. In the steady state conditions of catalysis, 
however, the surface of the alloy could be altered by the interplay of various processes, such 
as molecular adsorption and desorption and chemical reactions happening on the surface. 
Furthermore, in our thermodynamic equilibrium we only accounted for the presence of oxy- 
gen, while neglecting the reducing effects of ethylene. Our work can therefore be considered 
only as a first, albeit necessary, step toward understanding the full catalytic cycle. 
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